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Abstract. We show that in multidimensional gravity vector fields completely 
determine the structure and properties of singularity. It turns out that in the presence 
of a vector field the oscillatory regime exists in all spatial dimensions and for all 
homogeneous models. By analyzing the Hamiltonian equations we derive the Poincare 
return map associated to the Kasner indexes and fix the rules according to which the 
Kasner vectors rotate. In correspondence to a 4-dimensional space time, the oscillatory 
regime here constructed overlap the usual Bclinski-Khalatnikov-Liftshitz one. 
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1. Introduction 

The wide interest attracted by the homogeneous cosmological models of the Bianchi 
classification pP relies over all in the allowance for their anisotropic dynamics. Among 
them the types VIII and IX stand because of their chaotic evolution toward the initial 
singularity [2 . The stochastic properties outlined by such models were longly studied 
in order to describe their detailed features j3 IH El and in recent years a new line 
of research has been developed to establish if the chaos is covariant and survives in 
each system of space-time coordinates [3 El 13 EH HI] (see also references therein. 
For approaches via Ashtekar variables which, providing a non-negative Hamiltonian 
function, could have interesting applications in the analysis of chaotic anisotropic 
cosmologies, see the formulation presented in [TS]). 

The cosmological interest in the Bianchi types VIII and IX universes dynamics comes 
out because they correspond to the maximum degree of generality allowed by the 
homogeneity constraint; as a consequence it was shown (THl HH EEH1 EE! CZ| that the 
generic cosmological solution can be described properly, near the Big-Bang, in terms 
of the homogeneous chaotic dynamics as referred to each cosmological horizon (recent 
interest on a special homogeneous model like the Taub one arose because in ^S] was 
shown that it admits an accelerating stage in the dynamics as the actual Universe seems 
to do). 

It is just the dynamical decoupling of the space points which takes place asymptotically 
that allows such an inhomogeneous extension; indeed the increase of the spatial gradients 
is controlled by the oscillatory behavior of their time dependence, so that their effects 
are negligible at horizon size However the correspondence existing between the 

homogeneous dynamics and the generic inhomogeneous one holds only in four space-time 
dimensions. In fact a generic cosmological inhomogeneous model remains characterized 
by chaos near the Big-Bang up to a ten dimensional space-time [HI EUl EH 1221 EHj while 
the homogeneous models show a regular (chaos free) dynamics beyond four dimensions 

[2HEEI. 

Here we address an Hamiltonian point of view and show how the homogeneous models (of 
each type) single out, near the singularity, an oscillatory regime in correspondence to any 
number of dimensions, as soon as an electromagnetic field is included in the dynamics 
(for a study of the dynamics induced by an electromagnetic field in the framework of 
the quasi-isotropic solution, see [2E])- In fact the presence of an electromagnetic field 
restores a closed potential domain for the point-universe evolution; thus, even if the 
considered model does not posses a sufficient numbers of potential walls to determine 
the oscillatory regime, its asymptotic evolution is dominated by the matter field and 
randomizes within the corresponding billiard (for a connected analysis which studies 
the disappearance of chaos when a scalar field is included in the dynamics, though in 
presence of a p-form, see |2*7j). 

We stress how the electromagnetic billiard and its associated Poincare return map 
coincide with usual ones in four dimensions. However for higher dimensional spaces 
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they differ from the corresponding vacuum dynamics; in particular the map we derive 
outlines an interesting dependence on the number of dimensions. A relevant out coming 
of our analysis consists of fixing the law for the rotation of the n-independent oscillating 
directions. 

The procedure here followed calls attention to be extended to more general contexts 
like those studies on the appearance of chaos in superstring dynamics [23123- More 
precisely, in Section 2 we review the Hamiltonian formalism for the General Relativity; 
in Section 3 we derive the Kasner solution for homogeneous models. In Section 4 we 
introduce the Kasner parametrization in our dynamical scheme. Finally in Section 5 
we derive the solution in the Kasner approximation, and in Section 6 we show how this 
dynamical system can be described in terms of the return map for the Kasner indexes 
and the law for Kasner vectors rotation. In Section 7 brief Concluding remarks follow. 



2. Hamiltonian formulation 



We start by reviewing the dynamical framework of our analysis. 

Let us consider a vector field = (<p,A a ) ,(a = 1,2, ...,n) and adopt for the 
metric the standard ADM representation |37H I3*T] 

ds 2 = N 2 dt 2 - g af3 (dx a + N a dt) (dx? + N p di) . (1) 

Then the action, which describe the dynamics of the model, takes the form 



n Q/3 ^ Q/3 + 7r a -A a + <pD a ir a - NH - N a H a 



(2) 



where 



H = i= {nX - ~~J TO 2 + \dap^ + 9 (jFapF* - fl) } , (3) 

H a = -V^ + 7r^F a/3 , (4) 

denote respectively the super-Hamiltonian and super-momentum; here F a p = 
dpAa — d a A/3 is the electromagnetic tensor, g = det(g a p) is the determinant of the n- 
metric, R is the n-scalar of curvature constructed by the metric g a /3 and D a = d a + A a . 
n a and U alS are the conjugate momenta of the electromagnetic field and of the n-metric 
respectively ; they result to be a vector and a tensorial density respectively of weight 
1/2 because their expressions contain ^/g, and are defined via the relations: 



U a(3 = V_^(K Q ^ - g aP Tr{K)) (5) 



7T 



Ei (»4y* _ n> f A (6 ) 



(here K af3 denotes the extrinsic curvature in the synchronous frame) 



Benini, Kirillov and Montani 



4 



When varying this action with respect to the lapse function iV we obtain the Super- 
Hamiltonian constraint 

#o = 0, (7) 

while its variation with respect to ip provides the constraint d a n a = 0. Since 
we are dealing with an abelian vector field (i.e. corresponding to an abelian group 
of symmetry like an electromagnetic field does) whose sources (charged particles) are 
absent, it is enough to consider only the transverse (or Lorentz) components for A a and 
7T a . Therefore, we take the gauge conditions ip = and D a 7i a = and this will be 
enough to exclude the longitudinal parts of the vector field from the action. 
It is worth noting how, in the general case, i.e. either in presence of the sources, either 
in the case of non-abelian vector fields, this simplification can no longer take place in 
such explicit form and, therefore, we have to retain the term <p(d a + A a )n a in the action 
principle. 

In what follows we consider the behavior of homogeneous cosmological models in the 
asymptotic limit toward the initial singularity. 



3. Homogeneous cosmological models: basic equations 

When considering the homogeneity constraint, the whole spatial dependence of the 
models can be integrated out from the action and the dynamical variables become 
only time dependent. In general neglecting the total Hamiltonian the spatial 
derivatives contained in F a/3 and R, is equivalent to the so-called (generalized) Kasner 
approximation. 

When going over the homogeneous case, we choose the gauge N a = and, within 
the Kasner approximation, we get equations of motion for the vector field having the 
form 

d N 

|^ = 0. (9) 

The field equations which describe the n-metric dynamics read as follows 
d 2N f 1 



dt I 



{n a/3 - — [ g aP Il7 ( } , (10) 



^ = -T-^V (U) 



^rr - -JL 

This dynamical scheme is completed by adding to the above Hamiltonian equations 
the Super-Hamiltonian constraint (jJJ). 
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4. Kasner parameterization 

To develop the below analysis, it turns out very convenient to adopt the so-called Kasner 
parameterization which is based on the metric and conjugate momentum decomposition 
along spatial n-bein: 

where the n-bein is chosen in such a way that the matrix p ao = diag (pi, . . . ,p n ).(for 
a discussion of this hamiltonian structure in correspondence to an inhomogeneous 
multidimensional model see this diagonal form of the conjugate momenta is a 

consequence of requiring the canonical nature of the adopted transformation. We also 
define a dual basis L a a = g a H%, such that L a a l\ = 5 b a and L a a l% = 5%. 
To face our goal we have to project (fTUj) along the Kasner vectors defined in (fP2j) : in 
this way we get the following dynamical system: 

(4'«) = ^-^?4 (13) 

Ua^lbj + (^J^aJ =0, a^b. (14) 

Here (LJ b ) = L°l ba denotes the ordinary vector product, treating the vector 
components as Euclidean ones. 

In close analogy with above, from equation (fTTj) . we find the additional equations 

\ at J 2^/gpa-pb 
In particular, we see that (|TB|) already contains ()14jl. 

By combining together both such systems, we obtain ()15|) as the first independent 
equation; the equation for the Kasner vectors takes place in the form (for the sake of 
simplicity we neglect the vector index a) 



dt a ^ 



We want to put in evidence the oscillatory regime that the bein vectors possess; to 
better analyze this dynamics let us distinguish scale functions in the following way 

l a = exp (<f/2) £ a . (18) 

L a = exp (-q a /2) C a . (19) 



being £ a the so called Kasner vectors. 
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Thus instead of (|15J1 and (|17|) we rewrite 

^p a = -T^^exp (q a ) , A a = (ttX) , (20) 




|? a = ^ba-— T Eft). (2i) 

|4 = E ^ exp 4 = E ( 22 > 

^ 2 ^ ^ Pa - Pb v y ^ Afo (pa _ 

The last two equations are obtained by substituting variables (jl8j) and then 
projecting (fTTjl on the vectors ; when the index c is equal to a, we get ([2T]) . otherwise 
we get ()22|) (in which we neglected a term vanishing in the exact Kasner-like solution, 
and therefore, in our analysis, of higher order). 

To this system we should also add the Hamiltonian constraint 

ff0 .0^{ & *-^( &o ) 2 + i5>^}. (23) 

To require that the quantities A a are constants implies that vectors £ a do not depend 
on time in turn and therefore no Kasner vectors rotation takes place (see section 6); 
such a situation corresponds exactly to the n-dimensional Kasner behavior [201 EH EH]. 
Since the rotation of Kasner vectors is induced by an higher order term for g —>■ with 
respect to the pure Kasner dynamics, we can assume that A a are near to be constant and 
therefore near the singularity (j2*Uj) . (|2T|) give, in the asymptotic limit to the cosmological 
singularity, the billiards [HI (HUE] on {n — 1) —dimensional Lobachevsky space, exactly 
like in the 3-dimensional mixmaster case [33]. By other words, in our scheme the 
evolution is not simply Kasner like, but we get a dynamical picture in which the point 
Universe moves according to a piecewise Kasner solution. The features predicted by the 
two equations (f2~U|) and (f2~Tj) will be discussed in section 5 and 6. 

We remark that the validity of the Kasner approximation is based on the possibility 
to neglect the n-dimensional Ricci scalar with respect to the term containing time 
derivatives. In the Kasner solution this picture holds along the whole system evolution 
to the initial singularity [33] ; instead in the piecewise Kasner solution the n-dimensional 
Ricci tensor is negligible only for finite time intervals ending with a bounce against 
the potential terms arising from the spatial curvature [2] . However in both these cases 
the validity of the Kasner approximation is confirmed by a large number of theoretical 
and numerical investigations, concerning homogeneous and inhomogeneous cosmological 
models near the singularity [21 El EJ 13 EJ UHl HS1 [HI IIEl CZl EH 

Our scheme relies on the choice of a synchronous (or Gaussian) reference frame 
which corresponds to have N — 1 and iV a = 0. This system has many interesting 
properties, among which out stands its geodesic nature [33]; in fact the line x a = 
constant results to be geodesic of the manifold and since the normal to spatial 
hypersurfaces reads n M = (1,0) free particles are co-moving to this system. ^From 
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a cosmological point of view the Gaussian frame is relevant because the galaxies are 
(almost) free particles in the actual Universe and therefore our phenomenology relies on 
a synchronous gauge. However the results obtained in this work have to remain valid 
for other gauge choices in view of the general covariance. 



5. Kasner solution 



To analyze the time evolution of the Kasner vectors and to obtain their rotation law, it 
result to be convenient to project them on the time independent quantities n a and to 
decompose them in the two parts: 

L = 4n + L±; l A \ = = o. (24) 

Here, for the sake of simplicity, we use the vector notation which is useful when treating 
the components of these vectors as Euclidean ones. 

Hence we can split ()22|) into the two independent components 



dt b^a. h (Pa - Pb) 

|^ = E r7^ A V ~L = At (t) i b± . (26) 

dt b+a. h (Pa ~ Pb) 

— * 

Since above the matrix A% does not depend on £ a ±, then we get the formal solution 

4i (*) = T exp 1 1 A a b (t') dt' | £ b± (t ) , (27) 

The remaining equations (|25[1 together with (J2Up. ()2ip. and ()23p provide a self-consistent 
dynamical system. 

Since the Kasner solution corresponds to neglect the contributions of the n- 
dimensional Ricci scalar to the dynamics, then to get such a behavior we have to take 
the limit in which all the terms exp (q a ) become of higher order; under this assumption 
we get the following simplified dynamical system 

p a = const, X a = const, £ a ± = const, 

feh=f (Pa-^E^), (28) 

E^-^r(EPa) 2 + |Ee^ = o, 

whose solution, in the gauge N — 1 and toward the cosmological singularity (g — > 0), 
takes the form 

gap = E t 2Sa m > s a =l-(n-l)-^, (29) 

a l^bPb 

where the Kasner indexes s a satisfy the identities 

E^ = E^ = i- (so) 



Benini, Kirillov and Montani 



8 



Let's take the Kasner indexes in the increasing order 

si < s 2 < . . . < s n ; (31) 

in this way we always have si < and s n > > 0. Therefore the largest increasing 
term (as t — > t Sl — > oo) among the neglected ones comes from si and it is to be taken 
into account to construct the oscillatory regime toward the cosmological singularity. 



6. Billiard representation: the return map and the rotation of Kasner 
vectors 

To construct the oscillatory regime we retain just the leading term corresponding to 
exp (g 1 ), so the field equations (}20|) and (J2T|) rewrite 

d ~ 

^p a = 0, a ^ 1, (32) 
9 N T2 / i\ 

a 2at / i \ 

= — [Pa ~ l^Pb ■ 

dt ^/g\ n - 1 Y / 
The first of equations (f3~2"J) gives Ai = const, while the second admits the solution 

A a (Pa - Pi) = const. (33) 

The remaining part of the dynamical system allows us to determine the return map 
governing the replacements of Kasner epochs (i.e. intervals of time during which the 
evolution is Kasner like) 

S'l = , < = Sa+ \~ 2S \ (34) 

Ai = Ai, \' a = Xa(l- 2 ^T 1 ^ 1 ) , (35) 

\ (n-2)s a + nsi J 

Defining the quantities k a = p a — Pb we find for them the following iteration 

law 

k[ = -ki (36) 

K =k a + (37) 

E*'=E* + ^2*i (38) 
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Our analysis is completed by investigating the rotation of Kasner vectors £ a through 
the epochs replacements; by (}2l)j) we get the following system 



8r = JM*^, |f 1± = o, (39) 
9t Ai (p« -pi) dt 

admitting the integral 

Ai4i - A/ii. = const. (40) 
Putting together (j2U), (|35|) and fj39j) we arrive to the final iteration law 

4 = 4 + <T«/l, ^a = ^f^ = 

Ai 

= -2, (n 7 1)Sl ^. (41) 
(n - 2) s a + nsi Ai 

which completes our dynamical scheme. 

Thus the homogeneous Universes here discussed approaches the initial singularity 
being described by a metric tensor with oscillating scale factors and rotating Kasner 
vectors. Passing from one Kasner epoch to another one, the negative Kasner index si is 
exchanged between different directions (for istance l\ and £2) and at the same time this 
directions rotate in the space according to the law (j41j) . The presence of a vector field is 
crucial because, independently on the considered model, it induces a (dynamically [Iff] ) 
closed domain on the configuration space. The amplitude of ^"-oscillations increase 
approaching the initial singularity and their minimum value approaches —00 

In correspondence of this oscillation of the scale factor the Kasner vectors i a rotate 
and, at the lowest order in q a , the quantities a a remain constant along a Kasner epoch; in 
this sense the vanishing behavior of the determinant g approaching the singularity does 
not affect significantly the rotation law (|4ip. We conclude by stressing that duration of 
a Kasner epoch decrease as the singularity is approached in this scheme. 



7. Concluding remarks 

By the study above developed, we have shown how any multidimensional homogeneous 
cosmological model acquires, near the singularity, an oscillatory regime when an 
electromagnetic field is included in the dynamics. This result is a consequence of the 
capability that a vector field has to generate a billiard configuration in the asymptotic 
evolution; such a "billiard-ball" representation of the universe dynamics coincides with 
the BKL (Belinski-Khalatnikov-Liftshitz) piecewise approach only in the 4-dimensional 
space-time while in higher dimensions new features appear. In particular the obtained 
oscillatory regime characterizes all the homogeneous models, disregarding their potential 
term; furthermore the map by which the Kasner indexes evolve acquires a direct 
dependence on the dimensions number. 

A valuable issue of our Hamiltonian dynamics relies on fixing the rule of the Kasner 
vectors rotation. Such a law of rotation is relevant to connect the Cauchy problem with 
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later stages of the system dynamics. More precisely the rotation of the Kasner vectors 
is sensitive to the boundary conditions on the "matter" fields and has to be taken into 
account when using the studied homogeneous models within a cosmo logical "picture". 

Once the spatial gradients are taken into account, our result can be extended to a 
generic inhomogeneous cosmological model, in the same spirit as the Bianchi type VIII 
and IX model oscillatory regime is upgraded in four dimensions. 

Such an extension will reliably show how, in the presence of a vector field, the 
generic cosmological solution is described, in correspondence to any dimensions number, 
by an oscillatory approach to the Big-Bang. 

Investigations in this direction, as well as extended to the more general frameworks 
of superstrings [23 I2H] and brane dynamics [33] have to be subject for further 
developments. 
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